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[57] ABSTRACT 

Two alternative methods are disclosed for digital re- 
duction of synthetic aperture multipolarized radar data 
using scattering matrices, or using Stokes matrices, of 
four consecutive along-track pixels to produce “aver- 
aged” data for generating a synthetic polarization im- 
age. 


5 Claims, 4 Drawing Sheets 
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FIG. 2 
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Intensity Corresponding to a Four-Look Pixel 
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DATA VOLUME REDUCTION FOR IMAGING 
RADAR POLARIMETRY 

ORIGIN OF THE INVENTION 5 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provision of Public Law 96-517 (35 U.S.C. 
202) in which the contractor has elected not to retain io 
title. 

TECHNICAL FIELD 

This invention relates to a method for digital reduc- 
tion of synthetic aperture multipolarized radar data 15 
while still allowing full polarimetric utility of the data. 
The technique results in about 12.8-fold reduction in 
volume to allow multipolarized radar data to be used by 
a university and small industry users. ^ 

BACKGROUND ART 

Recent radar measurements show that different 
ground locations respond differently when the polariza- 
tion of either the receiving or the transmitting antenna 25 
is varied (see Dino Guili, “Polarization diversity in 
radars,” Proceedings of the IEEE, Vol. 74, No. 2, Feb- 
ruary 1986). In 1985, a polarimeter was flown on a 
NASA aircraft, recording data which were subse- 
quently processed and stored at the Jet Propulsion Lab- 30 
oratory. One difficulty inherent in this experiment was 
the large storage necessary for each of the images; only 
a limited number of data sets could be stored. Further- 
more, the computer time necessary to synthesize a pic- 
ture using an arbitrary transmit and receive polarization 35 
is mostly devoted to cumbersome data transfers. 
Clearly, the volume of data needs to be reduced in order 
to provide the user with a more flexible investigation 
tool, yet the data volume reduction must not impact the 
noise level by introducing additional error. Some theo- 
retical concepts will first be presented for clarity and 
review, as well as a description of the operations neces- 
sary to generate a synthetic polarization image from the 
original data sets. Then two new methods of data com- 45 
pression will be described. 

STATEMENT OF THE INVENTION 

In a first embodiment of the invention, the scattering 
matrices of four consecutive along-track pixels are “av- 50 
eraged,” while in a second embodiment, the Stokes 
matrices of four consecutive along-track pixels are “av- 
eraged.” In both cases, the average matrices are stored 
as a reduced data set. Each approach in terms of data 
volume reduction and in terms of errors introduced in 55 
the synthesized images substantially reduces the size of 
mercury required for each of the images, and the num- 
ber of operations is also reduced, thus reducing data 
processing time. This is accomplished in one approach 
by compressing scattering matrices, and in the other by, 
compressing Stokes matrices. The reduced data set is 
smaller in the first approach than that created by com- 
pressing Stokes matrices. However, greater error is 
introduced by compressing scattering matrices: these 65 
errors may be as great as 10% to be compared to typi- 
cally less than 10~ 3 when the compression algorithm 
operates on the Stokes matrices. 


2 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a flow chart of a prior art process 
for generating a synthetic polarization image from a set 
of imaging radar polarimeter data. 

FIG. 2 illustrates a flow chart of a first new method 
(scattering matrix approach) for generating a synthetic 
polarization image with data volume reduction. 

FIG. 3 illustrates a flow chart of a second new 
method (Stokes matrix approach) for generating a syn- 
thetic polarization image with data reduction. 

FIG. 4 is a typical histogram of the Stokes matrix 
elements F mn used as an example in determining the 
preferred method for optimizing the Stokes matrix 
method of data reduction illustrated in FIG. 3. 

DETAILED DESCRIPTION OF THE 
INVENTION 

The expression for the electric field of an electromag- 
netic wave propagating along the z axis is: 


E(z,t) = Re 


Eriiz.i) 

Ey(z,t) 


Re(h e'Xte-w)) 


(1) 


where, 


h = 


oh 

a v e — 


|exp(/B//), and 5 = 8// — 8 V 


where h is the polarization vector, Ejyand Eyare the 
horizontally and vertically polarized components of the 
electric field with respective magnitudes a// and a [/and 
arugments S h and 8 y. 

Another way to describe the same electric field is 
given by the Stokes parameters, defined as follows: 


Go=a[p- + av l 


G\—a}y — a\P- 


Gi = 2 a j]a, cos (6) (2) 

G3=2a#aj/sin (8) 


The Stokes vector is the vector whose components are 
the four Stokes parameters. If Go 2 = Gi 2 + G2 2 + G3 2 , the 
electromagnetic wave is said to be fully polarized. 

If the backscattering phenomena is assumed to be 
linear, isotrospic and homogeneous, the backscattered 
wave polarization vector h s may be expressed as 


h s = Sh( where 5 — 


S.xx S xy 
Syx Syy 


( 3 ) 


where hr is the polarization vector of the incident wave 
or transmitting antenna and S, the scattering matrix, is a 
2x2 complex matrix. If h r is the polarization vector of 
the receiving antenna, the complex amplitude of the 
received signal is: 


V=h r T Sh, (4) 

Thus, once the scattering matrix is known, a synthe- 
sized response may be computed for any desired config- 
uration of antenna polarization states defined by h r and 

hr. 
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Similarly, the 4x4 real Stokes matrix F relates the 
power of the signal received to the Stokes vectors defin- 
ing the polarization state of the receiving and transmit- 
ting antennas, G r and G ( . 


P— G r r FGt (5) 

If the reciprocatory principle applies, then the scat- 
tering matrix is symmetrical and so is the Stokes matrix. 
In that case, the Stokes matrix can be derived from the 
scattering matrix as follows (the general relation can be 
found in “Light scattering by small particules,” H. C. 
Van de Hulst, Dover publication): 

Jxxxx — SxxSxx* Jxxxy = Sxx^xy* 

Jyyyy = RyySyy* Jxyyy = SxySyy* 

Jxyxy — S X yS X y* Jxxyy = SxxSyy* 

where the asterisk indicates the conjugate of the signal 
indicated. 


Fn = 0.25 ( Jxxxx + Uxyxy + Jyyyy) F 23 = 0.5 Re (J XX xy - Jxyyy)*'® 

F\1 = 0.25 (Jxxxx - Jyyyy) ^24 = 0.5 Im (-Jxxxy + Jxyyy ) 

Fi 3 = 0.5 Re ( Jxxxy + Jxyyy ) ^33 = 0.5 Re ( Jxxyy + Jxyxy) 

Fi 4 = 0.5 Im (-Jxxxy - Jxyyy ) ^34 = 0.5 Im (-Jxxyy) 

F 44 = 0.5 Re (-Jxxyy + Jxyxy) 

There is enough information in either the Stokes 
matrix, F, or the scattering matrix, S, to synthesize any 
polarization configuration. The scattering matrix yields 
a synthesized signal whose power is identical to the 
power obtained from the Stokes matrix as long as the 
waves are fully polarized. 

A wave is fully polarized if it can be expressed as the 
superposition of a horizontally polarized part (HP) and 
a vertically polarized part (VP). The polarization vec- 
tor describes the electric field as a combination of a 
vertical part and a horizontal part. Therefore this repre- 
sentation is always referring to a fully polarized wave. 
Since the scattering matrix is the linear operator associ- 
ated with a polarization vector, the scattering matrix 
representation of a scatterer assumes that this scatterer 
cannot introduce any diffuse component (part of a wave 
which is not polarized) in the backscattered wave. The 
backscattered wave has to be fully polarized because it 
is written as a polarization vector. The Stokes vector 
representation, as opposed to the polarization vector 
representation, allows for a diffuse component which 
can be estimated by the difference 
Go 2 — (Gi 2 +G2 2 +G3 2 ) (see H. C. Van de Hulst, supra). 
Consequently, the Stokes matrix representation of a 
scattering surface can include this depolarization phe- 
nomena. 

Three methods will now be described with a com- 
plete set of data from an imaging radar polarimeter 
flown in 1985 on a NASA CV990 aircraft. Out of the 
subsequent data set, three images were chosen which 
feature a variety of targets for experimental testing of an 
original (prior art) process illustrated in FIG. 1 and two 
new processes, one using a scattering matrix approach 
illustrated in FIG. 2 and the other using a Stokes matrix 
approach illustrated in FIG. 3. The original process 
requires about 128 megabytes of synthetic aperture 
image data storage for one scene (frame). The complete 
data set is comprised of 4.2 million one-look pixels, 
1024x 4 points in the along-track direction and 1024 


4 

pixels in the range direction. For each pixel, four com- 
plex elements are stored per scattering matrix. A, B, C, 
and D in blocks 10, 11, 12, and 13 represent four consec- 
utive along-track pixels. Since a complex number is 8 
5 bytes long, the storage requirements are 128 megabytes. 
If the reciprocity principle is assumed, then the scatter- 
ing matrix should be symmetrical and only the three 
different elements of the scattering matrix need be 
stored. In fact, this property is assumed when calibrat- 
10 ing the data. 

In the original process, it was necessary to generate 
one synthetic polarization image from this set of data in 
the following manner. The polarization vectors of both 
the transmitting and the receiving antenna (not shown) 
15 of the experiment are chosen. The received signal is 
synthesized for each pixel using equation (4), as indi- 
cated by blocks 15, 16, 17 and 18. The resulting powers 
are then computed in blocks 19, 20, 21 and 22. The last 
step consists of averaging the computed power of 4 
20 consecutive along-track points in block 23. This opera- 
tion reduces statistical variations and makes the element 
of resolution approximately square. Each resulting av- 
eraged pixel now has an associated power and the image 
?5 is ready for display, as indicated in block 24, namely as 
intensity relative to one four-look pixel. 

This process of synthesizing an image of arbitrary 
polarization requires about twenty minutes on a VAX 
785 computer. Future users of the data will want to 
30 interact quickly with the picture making a faster process 
highly desirable. Reducing the size of the data set de- 
creases the processing time. The consequently reduced 
storage requirements also allow the storge of more data 
sets where only a few were previously possible. 

35 Two different data reduction methods have been 
developed for reduced storage requirements using dif- 
ferent approaches referred to hereinbefore as the scat- 
tering matrix approach and the Stokes matrix approach. 
In each case, a matrix corresponding to four along-track 
40 consecutive pixels A, B, C and D is stored as the result 
of an “averaging” operation done directly on the matri- 
ces, and not, as in the original process just described 
above, done on the synthesized signal power. The two 
new methods are detailed in the following sections, and 
45 a comparative study of the results is presented in a 
subsequent section. In each approach, the four complex 
elements A, B, C and D of the scattering matrix are 
stored as indicated at the top of FIGS. 2 and 3 in blocks 
10, 11, 12 and 13. 

50 THE SCATTERING MATRIX APPROACH 

An “average” of the scattering matrices is computed 
rather than an average of the synthesized signal powers: 
the “average” scattering matrix corresponding to one 
55 four look pixel is then the only information needed to 
synthesize a picture. FIG. 2 is a flow chart which de- 
scribes the process. For convenience in understanding 
this first embodiment and the second to be described 
afterwards, the same reference numerals ae employed 
60 for the same elements or blocks as used in FIG. 1. First 
an average scattering matrix is computed in block 25 for 
every set of four consecutive points in the azimuth 
direction. This matrix is stored in 8 bytes of computer 
memory represented by block 26. Eight megabytes of 
65 memory are then required to store the entire reduced 
data set for one image. 

Each time it is desired to generate a picture from a 
reduced data set 27, the first step is to choose the trans- 
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mit and received polarization vectors in block 14. The 
synthesized signal is then computed using equation (4) 
and its power from P=VV*, but in this case the re- 
duced data set is of averaged scattering matrices of 
four-look pixels, as indicated by blocks 28, 29 and 30, 3 
thus reducing the number of operations required to 
synthesize a complete image. The process is faster than 
before because less data is being manipulated and fewer 
operations are being performed. Data reduction can 
only be effectively achieved by carefully choosing the 10 
“averaging” method. The scattering method will now 
be described in more detail. 

Let A, B, C and D be the four scattering matrices to 
be “averaged,” corresponding to four consecutive 
points in the azimuth direction and let S be the resulting 
matrix. 


A xx4xy L = 

1 &x: •Sxy j 

dyx^yy | 

1 By X Byy j 

CxxCxy Ip = 

j CxxPxy 

CyxCyy | 

1 CyxCyy 
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in synthesizing a complete image, instead of the four in 
the original method illustrated in FIG. 1. 

THE STOKES MATRIX APPROACH 

FIG. 3 is a flow chart which describes the second 
new method for data volume reduction. Using equation 
(5), each scattering matrix A, B, C and D in blocks 10, 
11, 12 and 13 is transformed into its corresponding 
Stokes matrix Fi, F2, F3 and F4 in blocks 31, 32, 33 and 
34. One interesting property associated with the Stokes 
matrix is that it yields directly the power given the 
Stokes vectors of both the receiving and the transmit- 
ting antennas (equation 7). Therefore the two following 
processes are equivalent; adding the synthesized power 
of the signal scattered from 4 different areas or adding 
the 4 Stokes matrices characteristic of the four areas, 
and then computing the power of the resulting signal 
from this composite matrix. Let Fi, F2, F3, F4 be the 
Stokes matrices associated with four consecutive pixels 
and, let Grand G r be the Stokes vectors of the transmit- 
ting and receiving antennas, respectively. The powers 
P, corresponding to the phase matrix F, are given by: 

Pi=G R T FiG, ( 8 ) 


In the output image, each pixel is associated simply 
with an intensity, and the absolute phase information of 
the four received signals can be discarded, although the 
relative phases re yet required. Therefore, the absolute 
phase of the scattering matrix can be set arbitrarily. The 30 
power in each element of the “averaged” matrix is com- 
puted as follows: 

|S«|=sqrt(A xx . 

‘‘Vt X * + B XV B XX* XX* + D, tx D xx*) ^ 

I S X y 1 = sqrt(A, v A x ..* 4- B,.B X> * ■- GxyCxy* -r D x . 

jD X y*) 

| Sy X | = sqrt(A^ x Ay X * 4- B ;V B, X * 4- C. X C VX * 4- D, x . 

Dy*') 0) 40 

I ®Jtvl — Sqrt(A^A^* H- ^yyByy* 4" Cjy/Cjjy* -f- - 

D^yDjy,*) 

where “sqrt” signifies “the square root of ...” In the 
following, the phase of S** is set to zero. As discussed 45 
before, this is legitimate as long as the phases of the 
other terms relative to S xx are preserved. 

This method for calculating the phase puts more 
weight on the vectors of large amplitude which are 5Q 
presumably less sensitive to noise. The information to 
be coded consists of four amplitudes and three phases 
and can be stored in 8 bytes as follows: 

2 bytes for the maximum amplitude stored as a man- 
tissa and an exponent. 55 

3 bytes for the three remaining amplitudes. 

3 bytes for the three phases. 

After reduction, the data consists of 1024 X 1024 pix- 
els and 8 bytes per pixel, resulting in an overall reduc- 
tion ratio of 16. The reduced data set may then be pro- 60 
cessed to synthesize a complete image, using one “aver- 
aged” scattering matrix per pixel, as indicated in block 
28 and then proceeding in blocks 29 and 30 by comput- 
ing only one signal using equation (4) and then comput- 
ing power by the equation P=VV* to produce the 65 
same signal in block 24 as in the original method. Thus, 
in reducing data by averaging the four scattering matri- 
ces in block 25, there is only one signal processing path 


The four-look averaging process consists of adding 
the powers of four consecutive pixels as indicated in 
block 35. The resulting P can be expressed as follows: 

P^%P^G r T FiG t =G r T (XFDG t (9) 

Therefore, it is possible to add the four Stokes matri- 
ces corresponding to four consecutive points in the 
azimuth direction to form a single Stokes matrix, as 
indicated in FIG. 3 by block 36. This process is equiva- 
lent to a four-look averaging operation. The resulting 
Stokes matrix corresponding to a four-look pixel is then 
stored in block 37 in a compressed form as the reduced 
data set. 

To form an image from this reduced data set stored in 
block 37, the transmit and received Stokes vectors are 
computed (blocks 38) from the corresponding polariza- 
tion vectors 14 in equation (2). For each pixel, the 
Stokes matrix (block 39) is selected from the reduced 
data set 37 and the synthesized signal power is obtained 
in block 40 through equation (7). 

Each Stokes matrix in block 36 is a 4x4 symmetrical 
real matrix. It consists of 10 distinct elements, nine of 
which are independent (equation 5). The tenth, F22, can 
be computed fromm three of the nine elements as fol- 
lows: 

F-z2=F\\— /33— F44 (10) 

Fi 1 is the largest element and is always positive. It is 
coded in two bytes as a mantissa and an exponent. The 
eight other independent elements are scaled by F n and 
coded in 1 byte each. Therefore, the complete phase 
matrix can be stored using 10 bytes. The complete data 
set for one image (1024X 1024 pixels) is now stored in 
10x1024x1024 bytes. The overall reduction factor 
from the original data is 12.8. 

COMPARATIVE STUDY OF THE RESULTS OF 
THE DIFFERENT ALGORITHMS 

A testing procedure necessary to validate the differ- 
ent algorithms and to compare their performances com- 
pares a polarization signature obtained from the original 
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data set with the corresponding polarization signature 
computed from the reduced data set. A polarization 
signature is a surface in which height corresponds to the 
intensity of the signal received by an antenna transmit- 
ting a polarization defined by the x'and y coordinates of 5 
the point on the surface. The output of this program is: 


8 

corresponds to the exponent EXP and the second one 
corresponds to the mantissa MAN. 

EXP = E(ln(F\ \)/lnl), (12) 

MAN =F\ \/2 exp 


Error = \J 


J J U - l) 2 sin x d x dj, 
f f fism x d x d& 


where E(x) is the largest integer smaller than x. 

Since MAN is always less than 2 and greater or equal 
10 to 1, better precision can be achieved by storing MAN 
as: 


where I is the intensity of the pixel computed from the 
original data and I is the intensity computed from the 
reduced data set. The integration is computed over the 
set of all possible polarizations for the emitting antenna 
assuming the receiving antenna is the same as the emit- 
ting antenna (copolarized error) or assuming that the 
receiving and emitting antennas are orthogonal or 20 
crosspolarized (cross-polarized error). 

The test program outputs shown below indicate the 
difference between the original polarization surface and 
the one computed from the reduced data set for both 
reduction methods. 


ERROR (Eq. 11) 

Forest area 

Urban area 

Ocean area 


copol Scattering 

signa matrix 

7.19 E-2 

4.10 E-2 

8.80 E-2 


method 

Stokes 

matrix 

method 

4.06 E-4 

3.40 E-4 

3.14 E-4 

30 

35 

ERROR (Eq. 11) 

Forest area 

Urban area 

Ocean area 



cross- 

Scattering 

1.20 E-l 

7.20 E-2 

1.20 E-l 


pol 

matrix 





signa 

method 






Stokes 

4.18 E-4 

3.09 E-4 

4.25 E-4 

40 


matrix 






method 






byte(2)=£(254*(MAN-1.5)) (13) 

where byte(2) will range between — 128 and 127. EXP 
is stored directly as one byte: 

byte(l)=EXP (14) 

The eight remaining elements are first scaled by Fu 

Pij=Fij/F n (15) 

The absolute value of each P mn , where the subscripts 
mn designate a particular P y, is always less than 1. Two 
different strategies are possible for the storage of each 
term: storing P m „ directly or storing its square root. If 
P mn is digitized in 1 byte, the resul ting e rror is ±2~ 8 . If 
VP m „ is digitized, the error on VP,„„ is ±2 -8 , there- 
fore the resulting P mn is: 

(Vp^ ± 2- 8 ) 2 = P mn ± 2-7 

Comparing 2 -8 and 2~ 7 Vp i„„ shows that takin g the 
square root yields better results as long as VP m '„ is less 
than 0.5, i.e., Pmn <0.25. A typical histogram of each 
P m „ from actual data collected over the San Francisco 
Bay area is presented in FIG. 4. Pn, P 33 , P 44 are mostly 
larger than 0.25 while P 12 , P 13 , P 14 . P 23 , P 24 , P 34 tend to 
be smaller. 

In the case where the square root of P is stored, 


The Stokes matrix method yields smaller error than 
the scattering matrix method. The four-look averaging 
process on the original data set corrsponds to adding 
the power of four consecutive pixels and therefore is 
completely equivalent to adding the four corresponding 
Stokes matrices. The scattering matrices yield the com- 
plex field, not the power of the received signal and 50 
“averaging” four consecutive scattering matrices is not 
equivalent to a four-look average. The Stokes matrix 
representation allows for a diffuse component while the 
scattering matrix representation ignores it. For this 
reason, the Stokes matrix compression introduces 55 
smaller error in the data. The compression ratio is better 
in the case of the scattering method but the resulting 
errors are not acceptable. The Stokes matric algorithm 
was chosen and can be further improved as follows. 

OPTIMIZATION OF THE STOKES MATRIX 
METHOD 

The error computed above can be somewhat reduced 
by more carefully choosing the storage procedures. The 
method preferred is based on a probabilistic knowledge 65 
of the relative magnitudes of terms in the Stokes matrix. 

As discussed previously, Fu, the largest element of 
the Stokes matrix, is coded in two bytes. The first one 


45 w ( 16 > 

byte(0 = £(sign(P m „)* N \P mn \ *127) 

Where E(x) is the largest integer smaller than x and sign 
(x) = — 1 when x is negative and + 1 otherwise. When P 
is stored directly, 

by te(i) = £( 1 27 *P mn ) (17) 

Four different algorithms, each corresponding to the 
different numbers of square root elements were tested 
and the table below shows the differences between 
them. 



Algo. A 

Algo. B 

Algo. C 

Algo. D 

Pl2 

normal 

normal 

normal 

normal 

Pl3 

normal 

normal 

Square rt 

Square rt 

Pl4 

normal 

Square rt 

Square rt 

Square rt 

P23 

normal 

Square rt 

Square rt 

Square rt 

P24 

normal 

normal 

Square rt 

Square rt 

P33 

normal 

normal 

normal 

normal 

P34 

normal 

normal 

normal 

Square rt 

P44 

normal 

normal 

normal 

normal 
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For example, the square root of P 24 is stored for algo- 
rithm C and D according to equation (16). In algorithm 
A and B, P 2 4 is stored directly (equation (17)). 

The tests conducted on these algorithms consist of 
comparing the signature surfaces of a specific area, 5 
generated from the studied algorithm to the signature 
surfaces generated from the unreduced set of data. The 
test results, which are error percentages, are presented 
below. 


error 
(Eq. 11) 

Algo. A 

Algo. B 

Algo. C 

Algo. D 

forest 

CO 

4.06 E-4 

4.46 E-4 

2.80 E-4 

2.58 E-4 

area 

cross 

4.18 E-4 

3.71 E-4 

4.11 E-4 

3.48 E-4 

urban 

CO 

3.40 E-4 

3.95 E-4 

3.23 E-4 

3.25 E-4 

area 

cross 

3.09 E-4 

2.10 E-4 

2.13 E-4 

2.24 E-4 

ocean 

CO 

3.14 E-4 

3.07 E-4 

2.08 E-4 

2.09 E-4 

area 

cross 

4.25 E-4 

4.29 E-4 

2.51 E-4 

2.51 E-4 


The square root operation requires more processing /u 
when reducing the data as well as when generating a 
picture. A trade off has to be made between better pre- 
cision and a faster method. Algorithm C was finally 
chosen. The errors introduced by algorithm C are ^ 
smaller than when algorithm B is used. Algorithm D is L 
more complex and its results are comparable with those 
of algorithm C. 

The data reduction algorithm scheme finally adopted 
involves four square roots (Algorithm C): 30 

byte (3) = £(127 *Pi 2 ) (18) 

byte (4) = £(127*sign(Pi 3 )» N |Pi 3 | ) 

35 

byte (5) = £(127*sign(P 14 )» N \P H \ ) 
byte (6) = £(127*sign(/> 23 )* N |P 23 | ) 

40 

byte (7) = £(127*sign(P 2 4)* N |P 24 | ) 
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F(2,2)=F(1,1)-F(3,3)-F(4,4) 

CONCLUSION 

A method of reducing required storage space for 
multipolarimetric synthetic aperture radar data has 
been developed which preserves signal integrity. The 
data reduction operation reduces the required storage 
space by 12.8 and speeds the image synthesis process by 
a factor of 10 in time. The errors introduced in the 
output images are of the order of 10 -4 . Thus operations 
with polarimetric data are greatly facilitated and are 
within reach of even small research groups. 

. Although particular embodiments of the invention 
have been described and illustrated herein, it is reocg- 
nized that modifications and variations may readily 
occur to those skilled in the art. Consequently, it is 
intended that the claims be interpreted to cover such 
modifications and variations. 

We claim: 

1. A method for reduction of digital synthetic aper- 
ture multipolarized radar data while still allowing full 
polarimetric utility of said data, comprised of averaging 
and storing matrices corresponding to four along-track 
consecutive pixels, and using a set of average matrices 
thus stored to synthesize an image by computing a 
power signal for each average matrix corresponding to 
a four-look pixel after first choosing the transmit and 
received polarization vectors, thereby reducing the data 
set and number of operations required to synthesize an 
image for displays. 

2. A method as defined in claim 1 wherein the process 
of reducing said data set for one image by averging 
matrices corresponding to four along-track consecutive 
pixels comprises computing the average scattering ma- 
trix corresponding to each of said along-track consecu- 
tive pixels of an image, and wherein synthesis of an 
image from said reduced data set comprises first choos- 
ing transmit and receive polarization vectors, and then 
computing for each image pixel a synthesized signal 
from the equation 


byte (8) = £(127*F 33 ) 
byte (9) = £(127*P 34 ) 
byte (10) = £(127*P 35 ) 

RECONSTRUCTION OF THE STOKES MATRIX 

To reconstruct the Stokes matrix from the reduced 
data, the following operations are required: 

F(l,l)=(byte(2)/254+1.5)2 A r'e(l) 

F(l,2)=byte(3)*F(l,l)/127 
F(l,3)=sign(byte(4))*[byte(4)/127] I 2 *F(U) 
F(l,4)=sign(byte(5))*[byte(5)/127] 2 *F(1,1) 
F(2,3)=sign(byte(6))*[byte(6)/127] 2 *F(U) 

F(2,4) = sign(byte(7))*[byte(7)/ 1 27] 2 *F( 1 , 1) 


V=h r T Sh, 

where h r and h t are the desired receive and transmit 
polarization vectors, S is a two-by-two complex scatter- 
ing matrix and for each synthesized signal of a four-look 
pixel, computing its power P=VV* for each thereby 
producing a signal for each pixel having an intensity 
relative to a four-look pixel V. 

3. A method as defined in claim 2 wherein the four 
scattering matrices to be averaged are given by 

4 I ^ xx ^ x y Ijj I I 

I -'tyx^yy | | ByxByy j 


wherein the power in each element of the averaged 
matrix is computed as follows: 


c = 


60 


CyxGx 

CyxCvi 


u> = 


DxxD> cy 
DvxDvv 


F(3,3)=byte(8)*F(l,l)/127 
F(3,4) = byte(9)*F(l,l)/127 



IS^tl =sqrt(A xx . 

A*** + ByaBv.y* 4- C XX C XX * + Dy.yDa.y) 


F(4,4)=byte(10)*F(l,l)/127 


I | — sqrt(A^A A y* + 4- C + D r . 

y£>xy) 
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I Sy* | — Sqrt(Ay X Ay;t* + ByxByx* + CyxCyx* + Dyx- 
D y x *) 

\Syy\ — S^rt(A^Aj;j)* 4" YiyyQyy* *f CjyCjy* ~i* “ 

DjyDjy*) 

and the phase of S** is set to zero while the phases at 
other terms S XJ i, S^ x , and S w are preserved as follows: 

A.rg(S„)=0 

Arg(S X ji)= Arg(A x ^A xx * + B tJ ,B xx * + C X ^C XX * + Dv. 

^Dxt*) 


Arg(S^ x )=Arg(A^ x A XX* 4-B^jcB xx* 4-C_pxC XX- 

* + Dy X Dxx*) 

Arg(S yy) — Arg(Ayy\ XX * + B^>B XX * -t- CyyCxx- 

* + D w D.«*). 

4. A method as defined in claim 3 wherein each pixel 
has coded information said coded information includes 
four amplitudes, one for each of S xx , X yx and S yy , and 
three phases, one for each of S x ^, S^ x and S w in 8 bytes 
as follows: 

2 bytes for amplitude stored as mantissa and exponent 
for the first amplitude S xx , 
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3 bytes for the three remaining amplitudes S xv > S_ vx 
and S^y, 

3 bytes for the three phases, 
and the phase of S xx is set to zero so it need not be 
5 stored. 

5. A method as defined in claim 1 wherein the process 
of reducing said data set for one image by averaging 
matrices corresponding to four along-track consecutive 
pixels, comprises transforming each scattering matrix of 
10 four-look pixels into a Stokes matrix, and adding the 
four resulting Stokes matrices to produce a Stokes ma- 
trix corresponding to each four-look pixel, and wherein 
synthesis of an image from said reduced data set com- 
prises first choosing transmit and receive polarization 
15 vectors, transforming said chosen vectors into Stokes 
vectors of the transmitting and receiving vectors, and 
then computing for each time pixel a synthesized power 
signal from the equation 

20 A= G r T F,<;, 

to produce an intensity signal corresponding to one 
four-look pixel, where Grand G r are the Stokes transmit 
and receive vectors, P / corresponds to the Stokes matrix 
25 F,-, and the Stokes matrix F/is a four-by-four symmetri- 
cal real matrix. 
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